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We study the interplay effect of spin-orbit coupling(SOC) and optical lattice to the single-particle 
physics and superfluid-insulator transition in ultracold Fermi gases. We consider the type of SOC 
that has been realized in cold atoms experiments via two-photon Raman processes. Our analyses 
are based on the knowledge of full single-particle spectrum in lattices, without relying on any tight- 
binding approximation.We evaluate existing tight-binding models and point out their limitations in 
predicting the correct single-particle physics due to the missed high-band contributions. Moreover, 
we show that the Raman field (creating SOC) can induce band-gap closing in a two-dimensional opti¬ 
cal lattice, leading to the intriguing phenomenon of superffuidity-reentrance for interacting fermions 
at integer filling. We present the superfluid-insulator phase diagram in a wide parameter regime of 
chemical potentials and Raman fields. All these results are far beyond any tight-binding model can 
predict, and can be directly probed in current cold atoms experiments. 


As two typical potentials engineered in cold atoms, the 
spin-orbit coupling(SOC) and optical lattice both signif¬ 
icantly modify the single-particle dispersion and give rise 
to intriguing collective phenomena in interacting many- 
body systems[i|, Q. Their combination, i.e., spin-orbit 
coupled quantum gases in optical lattices, have recently 
attracted considerable attention in view of their exper¬ 
imental realizations through laser-assistant tunneling[^ 
1^, shaken optical lattices [3, Q, and two-photon Raman 
processes[§- 11|. A fascinating property of such system 
is that the effective flux in each plaquette is on the or¬ 
der of one flux quanta, large enough to reach the quan¬ 
tum Hall regime and enable the exploration of topolog¬ 
ical signatures [3-[l^. Theoretical studies have also re¬ 
vealed various interaction effects in such svstem [lB - E^ . 
and in particular, pointed out the possibilities of majo- 
rana fermions 2^, 221 and exotic spin textures fisl. [l^. 

To date, most studies on the spin-orbit coupled 
atomic gases in optical lattices have concentrated on 
the lowest-band 
approximations 
been rarely explored. Moreover, even for the lowest 
band(s), it is questionable whether the conventional 
Wannier wavefunction without SOC can still be used to 
construct the tight-binding model, as the SOC can in¬ 
duce coupling between different bands and the original 
Wannier basis could be problematic. With these motiva¬ 
tions, in this work we will go beyond the tight-binding 
approximation and lowest-band physics to explore the in¬ 
terplay effects of SOC and optical lattice. Our study will 
be based on the knowledge of full single-particle spec¬ 
trum, and thus will take into account all high-band con¬ 
tributions missed in previous studies. When come to the 
lowest band(s), this treatment also allows us to test the 
validity of tight-binding models used in literature 121-124 1. 

We consider the type of SOC that has been realized 
via two-photon Raman processes in optical lattices P- 111. 
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We exactly solve the single-particle spectrum, and point 
out that the existing tight-binding models have various 
limitations in predicting correct single-particle physics. 
This can be attributed to the SOC-induced excitations 
to higher bands, or equivalently, the improper usage of 
Wannier basis in the presence of SOC. In addition, we 
find a remarkable feature in the high-band physics, i.e., 
the Raman-induced gap closing in a two-dimensional(2D) 
optical lattice. With attractive interaction between two- 
species fermions, the gap closing leads to the reentrance 
of superfluidity from insulating phase at integer filling 
(two atoms per site). We identify the superfluid-insulator 
phase boundaries in a wide parameter regime of chem¬ 
ical potential and Raman field, and show rich density 
distributions in the trapped system as varying Raman 
fields. These results reveal the intriguing single-particle 
and many-body physics due to the interplay between 
SOC and optical lattice, which are far beyond any tight- 
binding model can predict. Our results can be directly 
probed in current cold atom experiments. 

Single-particle physics. We will first address the single¬ 
particle physics in the absence of interaction. The spin- 
orbit coupled atoms in a one-dimensional(lD) optical lat¬ 
tice can be described by the Hamiltonian (set h = 1) 

Hx = :^{Px - + Vo cos^ikLx), ( 1 ) 

2m 

here q is the momentum transferred through the Ra¬ 
man processes creating SOC, and LIr is the Raman field 
strength; Vq is the lattice depth; is the recoil momen¬ 
tum giving the lattice spacing a = tt/Zcl and recoil energy 
Er = k\l{2m). It is straightforward to exactly solve the 
eigen-system of Eq.(l) by writing the field operator as 
'0cr(x) = J2nk4‘nki^)'^nk, with the Bloch wave function 
Kki^) = -^J2G<k{G)e*V+G)^^ where n is the band 
index, k is the crystal momentum (s [—fcL,/cL]), and G 
is the reciprocal vector. In each {n, k} sector, the Hamil¬ 
tonian can be diagonalized in the {G, cr} subspaceji^. 
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In this work we consider the interesting case of fci = 2q 
, where an inhnitesimal Hr can generate a gap between 
the lowest two bands (see Fig.l). 

To capture the low-energy physics in deep lattices, 
a traditional way is to expand the field operator by 
the Wannier function of the lowest band, ipa-{x) = 
J2i^n=o,a-ix — Xi)cicr (* is lattice site), so to result in cer¬ 
tain tight-binding (lattice) model. For the present case, 
two types of such models have been employed [l2l - l2^ : 


{A) H = -tJ2{e -b -b/i.c.) 

i 

( 2 ) 

i 

{B) H = -t y^(c| ^c»+i,a -b h.c.) + Ofl y^(4tCi,; + h.c.) 


+it' 


( 3 ) 


Model (A) uses a shifted Wannier basis ujna{x) = 
(i/| = n = 0), i.e., the conventional 
Wannier function with a momentum shift for each spin; 
t is the hopping amplitude identical to that without 
SOCfl; / \uJo{x)\‘^ is the on-site spin-flip. 

This is equivalent to the idea of Peierls substitution under 
the gauge potential [1^ . In comparison, Model 

(B) uses the conventional Wannier basis ujn=o(x), with t 
the same as in (A) and t' = -^ f dxujo{x)-^u}o{x — a) 231. 
Model (A) and (B) are fundamentally different as they 
are based on distinct Wannier bases [25 1. 

In Fig.l, we compare the single-particle dispersion ex¬ 
actly solved from Eq.(l) (see Fig.l(al-a3)) with those 
from tight-binding models (A) and (B) ((bl-b3) and (cl- 
c3)) for several values of Hr. At Hr = 0, Fig.l(al) shows 
that the spectra of up- and down-spin are respectively 
shifted right- and left-side by a momentum q = fci/2, 
compared to that without SOC, and they cross each other 
at k = Q,zLkR. These shifts are due to the fact that PxtXz 
term in Eq. (1) can be gauged away by a unitary trans¬ 
formation U = With an infinitesimal Hr, energy 

gaps will open at these crossing places, see Fig.l(a2). 
These features can be well captured by model (A) us¬ 
ing the shifted Wannier basis, as shown by Fig.l(bl,b2). 
When increasing Hr further, however, model (A) fails to 
reproduce the correct spectrum. At large Hr, the spin 
tends to be polarized and the spectrum approaches to 
that of a spinless particle in a lattice potential, with the 
lowest two bands respectively featuring s- and p-orbitals 
(as indicated by Fig.I(a3)). On the contrary, model (A) 
always produces similar band structure as in small Hr 
case, merely with narrower band widths and larger band 
gap (Fig.I(b3)). For model (B), we see from Fig.I(cl- 
c3) that it cannot reproduce the correct single-particle 
physics for any value of Hr. At zero Hr, the momen¬ 
tum shift, as given by t' (^ ^[13), is hardly visible in 
Fig.l(cl); moreover, t' cannot be fully gauged away here. 
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FIG. 1: (Color online) Single-particle energy spectrum in a ID 
optical lattice with Vo = 5El, based on the exact Hamiltonian 
(first panel) and the tight-binding models (A) and (B) (sec¬ 
ond and third panels). The model plots have incorporated the 
constant shifts when reducing the exact Hamiltonian to the 
tight-binding ones[2^. The hrst, second, and third columns 
are respectively with Hr/El = 0, 0.2, 0.8. The color of the 
curves represents the weight of different spins in the corre¬ 
sponding wave function (red for d and blue for d). 


in contrast to the continuum case and model (A). Thus, 
model (B) is inapplicable when SOC is present. 

The breakdown of model (A) at finite Hr, as manifest 
by obvious discrepancies between Fig.I(a3) and (b3), is 
naturally related to the higher-band contributions missed 
in this model. This can be seen clearly in the large Hr 
limit, where the correct Wannier wavefunction to approx¬ 
imate the lowest band(s) should be with the spin trans¬ 
versely polarized and the orbit identical to the conven¬ 
tional uJo(x). Importantly, this wavefunction is a super¬ 
position of all levels of shifted Wannier bases {ojnaix)}, 
which means that many higher bands will contribute 
to the low-energy physics. Similarly, the breakdown of 
model (B) can also be analyzed. For instance, at zero Hr, 
the PxCTz term can induce excitations to higher bands. 
The resulted spectra turn out to be well described by 
the shifted Wannier basis, which is a superposition of all 
levels of conventional ones. 

We remark that the insufficiency or inapplicability of 
various tight-binding models is a direct consequence of 
the interplay between SOC and optical lattice. SOC 
couples the spin and orbital degrees of freedom, while 
in a lattice setup the spin flip is directly associated with 
the band excitation. If the SOC-induced band excita¬ 
tions are very serious, the tight-binding model based on 
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the original Wannier basis will no longer be valid. This 
mechanism is expected to hold for a general type of SOC 
and background setupf?^. 



FIG. 2: (Color online) Validity of tight-binding model (A), 
(a) rjv {v = 1, 2) as a function of Qr for given Vo = ?>El- The 
arrow locates the critical Raman field Q,% (see text), (b) Q.% 
and the corresponding ratio Q./t as a function of Vb. 


It is useful to identify a quantitative parameter regime 
for model (A) to be valid. Intuitively, one would com¬ 
pare with the width of lowest band(s) and the gap to 
higher bands. Nevertheless, a large gap is generally as¬ 
sociated with a narrow width, thus it is not immediately 
clear how the requirement of Hr changes with Vq. We 
thus carry out a numerical analysis to the band structure 
and then compare with model predictions. Based on the 
dispersion of model (A) [^, k = 0 and k = ±kL always 
have the same energy and each band has double maxima 
and double minima (Fig.l(b2,b3)). However, exact solu¬ 
tions show that as Hr increases, the energies at fc = 0 and 
k = become more and more deviated, and finally the 
double minima vanish in the first(lowest) band and dou¬ 
ble maxima vanish in the second band (Fig.l(a3)). To 
describe this feature, we introduce two quantities 


rh' 


\E,{k = 0) - E,{k = ±kL)\ 


(^ = 1,2), (4) 


here Wv is the band width from exact solutions. r]iy is 
always zero according to model (A), while is generally 
finite according to exact spectrum. 

In Fig.2a, we show that at given Vb, both rji and r ]2 
increase with Hr and finally saturate at unity. We de¬ 
termine a critical Raman field (H^) for the validity of 
model (A) by requiring both rji and 772 below 10%. In 
Fig.2b, we plot as a function of Vb and its correspond¬ 
ing {H/t)c-'We can see that for Vb G (5, 10)El, varies 
in a narrow region between 0.04 and 0.051?^, comparable 
to the lowest-band width (g (0.02, O.OO)^!!,) but much 
smaller than the gap to higher bands (g (2.43, 4.55)F1l). 
Remarkably, such weak SOC strength can already pro¬ 
duce considerable deviation of the band structure from 
the exact solutions. With such Vb and Hr, we find {H/t)c 
can range from 0.5 to 1.8, which sets the upper limit of 
H/t for model (A) to be valid. 

Band evolution in 2D. The evolution of ID spectrum 
as tuning Hr (shown in Fig.l(al-cl)) can induce even in- 





FIG. 3: (Color online) Band-gap evolution in 2D optical lat¬ 
tices. The energy dispersion can be decomposed as i?k = 
Ex{kx) + Ey{ky). (al,a2) ifa, and Ey for the lowest two bands 
at Vb = 5El and Hr = O.SEr. The band index is marked 
accordingly (see text). Red dashed lines show the prediction 
of model (A) to the band lx, 2x and ly. (b) The first (be¬ 
tween lx + ly and 2x -|- ly) and second (between 2x -I- ly and 
lx -I- 2y) band gap as functions of Hr at Vo = 5El. Dashed 
line shows the prediction of model (A) to the first band gap. 


teresting band evolution in higher-D optical lattices, such 
as in 2D lattices when adding another lattice potential 
along y—direction, say, V{y) = Vb cos^(fcL?/). Generally, 
the lowest band in x-direction and the lowest band in 
y-direction, noted as lx -I- ly, comprise the lowest band 
in 2D; using similar notations, 2x -V ly and lx + 2y re¬ 
spectively comprise the second and third lowest bands, 
see Fig.3(al)(a2). As increasing Hr, the gap between la; 
and 2a; increases, leading to non-trivial band-gap evolu¬ 
tions. As shown in Fig.3b, the first band gap opens at 
Hr = 0.175Er when the maximum energy of the first 
band (at (kx,ky) = (±A;i,±fci)) matches the minimum 
of the second band (at (iA;^, 0)), while the second band 
gap closes at Hr = 1.35£1l when the second and the 
third band touch at (0, EkL). Note that the tight-binding 
model (A) tend to overestimate the first gap as increasing 
H (see Fig.3b), and meanwhile, cannot predict the clos¬ 
ing of the second band gap. Particularly, the band-gap 
closing reflects the unique high-band physics for the spin- 
orbit coupled atoms in high-D lattices, which turns out 
to induce intriguing many-body physics as shown below. 

Super fluid-Insulator transition. We investigate the su¬ 
perfluidity of fermions loaded into a 2D optical lattice. 

The interacting Hamiltonian is 

driljl{r)Ho{r)ilj^,{r)-\-g j dr-ilj].{r)tjjl{r)tjji{r)tjj^{r)-, 

era' 

(5) 

Here Ho{r) = Hy, + Hy, with r = {x,y) and Hy = 
Py/(2m) -I- Vocos^(kRy); the bare interaction g can be 
related to the 2D binding energy Eh by the renormaliza- 
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tion equation g~^ = —+ Eb)~^ {S is the 
quantization area). 

Following the treatment of fermion superfluidity in op¬ 
tical lattices [28l - l^ . we consider the dominated intra¬ 
band pairing with the opposite crystal momenta and the 
same energy, which gives the pairing amplitudes 

^ A//^(-!/'n-k'0nk); A„k = ^ AgM^*. 

nk G 

( 6 ) 

with M^y. = EQ«i,-k(-Q)«Ik(Q + G). Here V'nk and 
®nk analogs of and as specified in 

the ID case. Given Eq.®, the thermodynamic potential 
1C = {H — fiN) can be calculated as 


fc = ^2 (-E-nk — fj- — \/(Enk — gY -f IA„k — An_kp) 

l[\,kx>0,ky 


S 

9 




(7) 


where Enk is the single-particle eigen-energy and /i is the 
chemical potential. The ground state can be obtained by 
minimizing /C with respect to {Aq}. In practice we have 
kept the lowest three G = (0,0), (2fci,0), (0,2fci,). We 
have checked that the results obtained will not be visibly 
altered by including more Ag into the minimization. 

Near the Superfluid(SF)-Insulator(IN) transition, all 
{Ag} continuously evolve to zero, thus one can expand 
K, as 1C= X^GG' AqCgg'Ag', with the matrix 


Cgg' =-<Igg' 

9 


- E 

n.kx >0,ky 


(^nk - ^n%)(^nk " ^^^k) 

2|Enk — mI 


( 8 ) 

The SF-IN transition can thus be determined by setting 
the determinant of G—matrix to be zero. 

In Fig.4a, we present the SF-IN phase diagram in 
the g — CIr parameter plane with given Vq = 5El and 
Eh = 3 El. a remarkable feature is that the IN phase 
with filling p = 2 appears as an isolated area surrounded 
by the SF bath in the phase diagram. This suggests that 
a system with fixed p = 2 will undergo a sequence of 
transitions from SF to IN and to SF again when increas¬ 
ing H/j, as explicitly shown by the density contour plot in 
Fig.4a. This can be attributed to two competitive effects 
induced by H/j: the spin-polarization does not favor the 
s-wave pairing and thus gives rise to the first transition, 
while the vanishing band gap (as shown in Fig.3b) in¬ 
duces pairing around the Fermi surface and gives the sec¬ 
ond transition. Here, we can see that the single-particle 
physics manifest itself well in the SF-reentrance of inter¬ 
acting many-body systems. 

Fig.4a also shows the SF to IN transition at half-filling 
(p = 1) as changing CIr. This is due to its two coop¬ 
erative effects: the spin polarization and the increasing 
(first) band gap, both of which disfavor the pairing super¬ 
fluid. Though the model (A) can also capture these two 



FIG. 4: (Color online) (a)Superfluid-insulator phase diagram 
with Vo = 5El and Eb — 3 El- Contour-density plots are 
shown for filling factors p — 2.2, 2,1.8, 1.5, 1, 0.5 (top to bot¬ 
tom). (b) Density profiles of atoms in a trapped system at 
CIr — 0.1, 0.7, l.di^L. Vo and Eb are the same as in Fig.4(a). 
We consider the trapping frequency u) = (27r)140HZ and the 
lattice spacing a = 425nm, giving Ro = •j2Ehj(vnoP') = 
lOO/fci. The total particle number is fixed as 5 x 10^. 


effects, we do not expect it give the correct SF-IN phase 
boundary at CIr > Er, as the band gap is considerably 
overestimated in that regime (see Fig.3b). 

The (p, CIr) diagram implies very rich density profiles 
in a trapped system as varying CIr. In Fig.4b, we show 
several typical density profiles for trapped atoms, ob¬ 
tained by applying the local density approximation with 
p{r) = ^(0) — V (r) at position r, and V (r) = mw^r^/2 is 
the external harmonic trap. We can see that as increas¬ 
ing CIr, the p = 2 and p = 1 IN phases (with flat-top 
density) respectively emerge near the center and edge of 
the trap. Increasing VLr further, the p = 2 IN phase gives 
way to SF and the flat-top structure disappears. This is 
directly related to the gap-closing and SF-reentrance dis¬ 
cussed above. Eventually, at large VLr the p = I IN phase 
will occupy a large region in the trap. 

Summary. In summary, we have studied the single¬ 
particle physics and the superfluid-insulator transition 
due to the interplay of SOC and optical lattices. We point 
out various limitations of existing tight-binding models in 
predicting the correct single-particle physics. We also re¬ 
veal interesting high-band physics that have been missed 
in previous studies, including the SOC-driven band-gap 
evolution and the resulted intriguing SF-IN transitions, 
which are far beyond any tight-binding model can pre- 
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diet. These results can be directly explored in cur¬ 
rent cold atom experiments. For instance, the single¬ 
particle spectrum can be probed by the momentum- 
resolved Bragg spectroscopy[3l|; the SF-IN transition 
and the reentrance physics can be detected by measur¬ 
ing the momentum space condensation fraction similar 
to the continuum case(^. or can be inferred from the 
density profile in a trapped system. 
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Supplementary Materials 


In this supplementary material we provide more details about how to exactly solve the single-particle Hamiltonian 
and to derive the two tight-binding models. 


I. Exactly solving the single-particle Hamiltonian 

We first write down the second-quantized Hamiltonian for (Eq.(l) in the main text): 

H = J dx 


( 1 ) 


(T=td 


We can expand the field operator in terms of Bloch wave functions of H^, 'ipaix) = J2nk4>nki^)'^n'k, with k lying 
within the first Brillouin zone (BZ) and n the band index. According to the Bloch theorem, one can further expand 
the Bloch wave function in terms of plane waves, 4>nki^) — Sg with G the reciprocal vector 

(integer of 2fci). Then the Hamiltonian H can be written as 


n,k^BZ G' G 

+ E E°ifc(^')E°'tfc(^){[^(^ + ^ + g)^ + 2^o]SgG' + -^SG±2kL,G'}'fpik'^nk 

nk&BZ G' G 

n,k£BZ G' G G' G 

In each {n,k} sector, H^k = T.gG' T.aa' <k(.G)ankiG')Hg^G'^ the is 

0 ... 0 


( 2 ) 


( ^\-N2kL hb/'t 

K)/4 4-(^-i)2fc. hi,/4 ... 


rro-o- _ 

^GG' — 


0 

0 


V 0 


0 

0 

^R 

0 


'-k+N2kL 


0 

0 

Hff 


^R 

0 

0 

=-l- 

-'k—N2kL 

Eo/4 

0 


0 

^R 

0 

VolA 

A 

-k-{N-l)2kL 


0 


0 

0 

^R 

0 

0 


'-k+N2kL 


(3) 


/ 


Here ef = ^ (k^g)^ + ^, and ^ (4-^-b2) x (4A^ + 2) matrix, where A^ is a cutoff number and in our practical 

calculation we take N = 50. A larger N will not visibly change the spectra for the bands that we are interested in. 


II. Deriving the tight-binding models (A) and (B) 

To derive the tight-binding model, we expand the field operator in terms of Wannier function of the lowest band, 
ipaix) = J2i^'n=o,tT{x — Xi)ciaj here i is the index of lattice site. Given the expression of the Wannier function 
uJn=o,a-ix — Xi), the second-quantized form of the Hamiltonian can be reduced to a certain tight-banding model. 


Model fAJ 


In this model, we take the shifted Wannier basis 

ojn^rix - Xi) = - Xi), 


(4) 




2 


with ^'1 = ±1 and n = 0. The parameters in the model (A) include: 
(1) the nearest-neighbor hopping term: 


= - J dxujQ^{x - + K) COS^ {kLx)\ujoa{x - Xi+i) 

J dxujo{x){-^ + Vo cos^ {kox))uio{x - a) 


= —e 


(5) 


^-hVbcos^ 
above equation, we have used the relation: 


here t = — f dxwo(x)(-^ + Vq cos'^{ k]^x))wo{x — a) is the hopping amplitude identical to that without SOC. To obtain 


1 + K) COS^ ^+Vo cos\kLx). 

2m 2m 


( 6 ) 


Therefore the phase shift in front of the Wannier basis is to gauge away the Px<Xz term, so that the resulted topping 
only differs by a phase factor compared to the conventional one (without SOC). 

(2) the on-site spin-flip term: 


= j dxijjl„{x)VLiiijJos{x) = VLii jdxe 


= J (ia;cos(2ga:)|wo(a:)P 

= n 


(7) 


In previous studies, the on-site spin-flip Q, is generally set to be identical to Hr, which is a good approximation only 
for deep enough lattice with very localized Wannier wave function. 

In addition, in writing the model (A) (Eq.2 in the main text), we have omitted the effective chemical potential term 
with 


ef = y dxujg^{x - + CqCOS^ {kLx)]uJo„{x - x,) 

f P^ 

= dxujo{x){-^ + VoCOS^{kox))ujo{x) 


( 8 ) 


In comparing the spectrum between exact solution and model solutions (Fig.I), we have incorporated above energy 
shift in Fig.l(bl-b3). 

For the lattice model (A), its dispersion can be obtained straightforwardly by using the Fourier transform Ci^ = 
-= J2k Finally, we get the spectrum of the lowest two bands in k-space 


Ei 2ik) = ± sin^(7rfc/fci) -|- 12^. 


(9) 


Model (B) 


In this model, we take the conventional Wannier basis which is spin-independent: 

iOna{x - Xi) = UJn{x - Xi), (10) 

with n = 0. The parameters in the model (B) include: 

2 

(1) nearest-neighbor hopping term induced by ^ -I- Vb cos^ (k^x), i.e., the first term of model (B) with amplitude 
t. 


(2) nearest-neighbor hopping term induced by 

iqi^a 


t. 




m 

m 


dxujQ^x - Xi)—uJo(x - Xi+i) 

d 

dxLOo{x)—uJo(x — a) 


= it'l'a] 


( 11 ) 
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where t' is given in the main text. For the deep lattice Vq ^ El, we find 



which is much smaller than t (see Ref.l cited in the main text) : 



This results in an invisible momentum shift for the up and down-spin spectra in Fig.l(cl). This analysis also holds 
for other values of q that is typically of the order of in realistic experiments (Ref.9,10,11 in the main text). 

(3) the on-site spin-flip term: 




= 


(14) 


Similar to model (A), in writing the model (B) (Eq.3 in the main text), we have omitted the effective chemical 
potential term with 


= / dxu}Q{x - Xi)[—pl +Vo cos^ (kLx) +—]u}o{x - Xi) 


= e^ + ^ 
* 2m 


(15) 


In Fig.l(cl-c3), we have incorporated above energy shift. 

For the lattice model (B), its dispersion can also be obtained straightforwardly by using the Fourier transform 
Ci„ = and finally we get the spectrum of the lowest two bands: 


Ei 2 ik) = ef — 2t cos(7rfc/fcL) ± 41'^ sin^(7r/c/fcL) -|- 


(16) 



